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1. » INTRODUCTION AND OVERVIEW 

Throughout most of your studies in -calculus, you 
will be concerned wrth the evaluation of . inte-grals . 
Various techniques suclws iiftegration by parts are 
^developed so that the ant iderivati ve can be t^termined 
and used to calculate the precise value o"f an integral 
accordijig to the fundamental theorem of calculus. 

However, there are some integrals which can not 
be evaluated by deterfninin^ an ant ider ivat ive'. The 
integrltl ' ' t 

re d:!f * 

which aris^ in probability theory is one of many 
examples: * For su^h -integrals , we must use 'other tech- 
niques, such as "numerical integration" to find approx- 
imate numerical values ^or such integrals. ' * 

^ The vary definition of the integral as the limit 
0?->j^^iemann ^um provides one numerical method for - 
cp^litTing a sum .jyhich approximates . 

^ (^f(x-J dx • / * ^ 

, where f(x) is a bounded function in the interval (a,b]. 
Other methods are the trapezoidal rule and Simpson's 
»rule. X, ^ * ' 

Th^.-^pQ^pose of this unit is to introduce these 
methods of ^numecical integration and des-e^^ibe- tRe 



related programming techniques required for h-a^vmg^* 
^ a computer carry out the , calculations . This uait 

requires very* little previous background or knowledge, 
about computers and their use. - * - ^ 

^Throughout the module, various references are 
cited and a comf>leie bibliography follows the main,^" ' 
text .* 'Students are encouraged to pursue furth^er 
readings that are suggested if they desire a deeper* 
♦ understaij^ding of "thii maternal. Most importantly, 
^ stucfents *are expected to'work eaqh of the several exer- 
ti,ses; for, these are" esTential, to, a geniune under- 
'Standing of the ideas developed. 

^ 2. APPROXIMATING INTEGRALS .USING RLEMANN SUMS 
X * - ^ ' , 

2.1 The Le^-Rectanj5;le ^Method ^ * * ' 

^ Our first attempt to solve the' problem of evaluating 

^ definite integrals involves using Riem^nn sums as 

approximations\. We consider a function f defined on a 
closed finite interval a £ x « b\nd a partition 
a = x^'<- <''x2 < ... < x^^_^ < = x,^-= b of thM^ 
interval. Fpr each j = 1 , 2 , . . . , N , let c be any 
point in the j^^^ subinterval of the partition, 

x.;^' < c. <.x.. 

The Riemann sum • ^ | 



N \ 



0^ 



J^f(cj).(x^-Xj.p = f(cj)-(x^-x 
^ f(c.2)-(x2-xp 



itegrai /^fCx) dx. 
j*-" term of this sum, f (c ■ ) • (x . -x. ' )• is described 

^ J J / 



is 9n approximation of the integral / f(x) dx . The 



geometrically in Figure' 1 as the area of the shaded 
recxa*ngle. The approximation of the integral that* is 
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Figure 1 . ^ 



determi,ned by' the partition is the sum^^of the areas, of 
N rectangles of heights f(Cj)— this area actually is 
negative when f*(c:) .is negative- -and widtlij^ j '^j 1-^*' 

In this general formula, the value of c- in the 
•subinterval Uj,j,Xj] can be choseri'arbitrarily and in 
fact so can the x^'s in the partition. T9 implement 
this (Tiethod on a computer, we need a systematic way of 
making these* -choices . One method is to choose a part.i 
tion in which' each -subinterval has the same width as 
each:6tiier subinterva\. This is ea^ijy accomplished 
by subdividing the interval (a,b] whose width is (b-a) 
into N equal' subintervals each with width (b-a)/N. A 
convenient and standard notation for the quantj^ty 
(b-a)/N is Ax. The. parti tion is then 

a = Xq, x^- = a + Ax, x^ f a + 2Ax, . . , 



so we have N subintervals 
i ' 

(a, a + Ax] , (a + Ax, a + 2Ax] , 
. . . , (a + (N-i)Ax-, b] - 

each of which has length Ax . ' 

Now if we choose c^ as the left 'endpoint of. each 
of these subintervals so that 



a , c. 



a + Ax , 



(N-l)Ax 



then we can form a Riemann sujn approximating the inte 
^ gral called th'e lef t-reatangl^ approxiwation. It is 
given by, the following expression: 

f(a)Ax + f(a + Ax)Ax ' 

+ f(a + 2Ax)Ax + f(a + (N-l)Ax)Ax. 

, This can be simplified to 

♦ ^ 

fl) l^i^) + f(a + Ax) + ... + f(a + (N-l)Ax))Ax 

This approximation is shown* in Figure 2. , | 

f(al' 




f(a+ (N-1)' 

' Area of rectangle I is f (a) • x 
Area of rectangle 1 1 is f(a+Ax)«Ax 
Area pf rectangle II I is f(a+2-Ax) -Ax 



Area of rectangle N is f (a+ (N-l)Ax); Ax 

Sum- of above is Formula (1) 



Figure 2. 



Let us uso the lef t -rectai^gl^ method to approxi- 
mate the integral • ' f 



!' 1 



I -m^ dX, 



^0 1+x' 

^ which in fact ye know to be 



arctah (1) - arctan (0) = j = \78559. 

Here^f(x) = 1/(1 + x^) . ' * 

Suppose we use N =^4.' Since b = 1 an4 a = t) we 
have Ax = 1/4^ Thus the partition, is • - 



a = 0, = 1/4, • . • 

^ x^ = 1/2, Xj = .3/4, b '= 1 . 

and the approximation (1) is ' 

tf(0) ^ hl/4) V £(172) + f(3/4)j / 1/4 ^ 

- (1 + 16/17 + 4/5 + 1V 253 1/4 , 

= (1 + .941 + . 800 + .'540).- iy4 

^ = .84 5. ' \ - \ 

Although we iatend to use, the computer jeo find better 
approximations by increasing the value of N, it is 
'important that yqv unders'tand hoW* t<j'tise these formulas 
with han^calculat ions also^. . * ^ > , 

, Now we are'read'y to devise a simple program" {which 
can easily be translated anto a languag^ like BASIC) 
that will compute an approximation to an integral by\ 
the left-rectangle inethod. As Forrtrtila (1) indicates', 
we ne^y-to have x range through the values. of a , a + Ax 
a +^2Ax, a ^ (N-l)Ay, i.e., x ranges from a to 

a + (N-l^Ax in steps bf Ax. We ^eed to calculate f(xy'^ 
for each^such x, and add up these terms to form a sum ^ 
S which we" then mOltipl/ by Ax to Compute the^ left- 
rectangle appi*oximatibn. .Suppose we use the computer^ 
variable IT (for "delta") to represent Ax; the following 



10 



program is probably what y6u would have designed. In 
this example program, we use the function f(x) = l/(x^+l 
which we define in the first line of code ^ 



Example Program 1 

100 DEF^NA (X) =.l/(xt2 +1) ' . 

110 *IN^UT A, B, Nl * • V 

120 LET D = {B-A](W ' 

160 LET S = 0,: * " 

140 FOR X =* A to (A +• (N-'r)*D) STEP a ' ^ , 

160 ^LET S = S + FNA (X)' 

* 160 NEXT X' . ' 

170 PRINT S*D; "IS THE LEFT-RECTAJJGLE APPRCJX." 

180 £ND ' ^ * > ' 

Note that^when the INPOT A,. B, N, ins truction 'is . * 
exectued, the computer will request 3 numerical valuas \ 
fOr^Ay B, and N respectively. '(In an interactive BASIC 
session, these must 'be separated, by commas when typed in. 
Line 170 computes S-D and prints this v^lue out along ' 
with th.e descriptive phrase.oenciosed in cfuotes^. ' ^ 

If we increase the, number of subintervaPs ^by using 
larger values of N, we expect to obtain a betteY approxi- 
mation to, the integral since the Ax will be smaller and 
the rectangle more closely matched to'-the .function". Ihis 
is so because the left-rectangle approximat^n being, a 
Riemann sum, tends to the vaiue o'f the integral as Ax ^s 
made smaller. ' ^ . , 

N6w, in Example Program 1, we used f(x) 5.1/(x'^ + l) 
which^as antiderivative arctan (x) . Thus, in .this- case 
can determine the'Cixact value of the integral 

* arctan (B) - arctan \(A) ' * 

and compare it with our numerical approximation. 

We cpn make such a comparison anytime we worTc with 
^an f(x) fdV wrf^clf we know the antiderivative. For ' * 



IX 



instance,. we could do this with ftx) — 4x^ - 3x" + 2x - 7 

whose antiderivative is x"* - + x^ - 7x,- We -can easily 

modify the code of- our E^cample Program 1 to handle this 
• new function >y changing tine 100 to 

100 DBF FNA (x) = 4*xt3 . 3*xt2 + 2*x - 7 , * 

To make it easy to compare our numerical approximation 
with the exact value we can add a line or two of code to 
calculate the exact value ^from the antideri^^ative and then 
print out that value. . ^fe could do this by adding the line 

175, PRINT Bt4 - Bt5 + Bt2 - 7*& - At4 +.At3 - At2 + 7*A 
e (Note this is not the most efficient way to*evaluate a 
polynomial, buf it is straightforu'ard and easy to write. 
Horner»s method would be much better; see UMAP Unit 263 ; 
Horner's Scheme and Related Algorithms , by Werner C, * ' 
Rheiftboldt. If we are going to perform such. an evaluation 
many time$, we might consider such improvements.) * ; 

Now that we have the program for the ^eft- rectangle 
method, we can perform many experiments.* For instance, 
we- can change the function we want to integrate. We can 
also vary the N (and thus .the__Ax), tp^examine _tLe J.'goodn£jsV* 
of our method. ' 



Exercise' 1 . 

^ ..R^n Example Program 1. for fCx> = 4x'^ . 3x^ + 2x - 7 for fixed 
valued of A and B but make N increasingly larger. -Copy ^ome of 
your results on notepapei' including the value^ of A, B, "N for ^ 
several -choices of N. Does the approximation gef closer to the 
«xact value of the definitS integral? 

: -■ i 2 ^ ^ 

^ Yw^migSf " have noticed" that if is only necessary to 
calculate once . fhe exact value of tjie ittcgrai . for given 
va'iufes' of N. .Thus line 175 really should not be-executed 
every time you run- the program f'or a new- N. If'^yaii pl«>ri,, ' 
tq^run many experiments, you should consider placing line " ' 
17S outside the "loop" qf your changes to N, for instance. 
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by doing the exact evaluation before you begin any approx- 
imations. This would also help you to determine if your 
program is' running correctly. This i s ^a standard tech- 
'nique*; nameJy^ testing a, program with l^wn rebults. 

2.-2 The Right-Rectangle Method ^ 

'Suppose instead of choosing c^ as tlie 1 eft - endpomt 
of each submterval, we choose as the r ight-endpomt 
of each submterval so that c j = a + Ax , c ^ = a + "ZAx, 
. , . , Cj^ = a + NAx = b. Then our approximating Riemann 
sum becomes^ . _ ' ^. 

(2) \ (f(a + Ax) +'f(a + 2Ax) + ... + f(b)]- Ax. 

As this formula indicates, we need to have x run through 
the values a + Ax, a +'2Ax, B, Calculat e f (xj * for 

each such x, and add' up these terms to form a sum S 
*\^^hich is theri'mulitplied by Ax to compute what is called 
the right-^rectangle approximation. 



Exercise 2 . 

^ Modify Example* Program 1 so that the sum S is calculated as the 
right-rectangle approximation. Type in this modified program on. a 
.xomputei:-terminal-with instruet"iorr~175 which prints out the~&xact ~ 
value of the f(x) dx for f(x) = 4x^ - 3x^ + 2x - 7," Run this 
for fixed values of A, B but use several increasingly larger values 
of N. Compare your results with those of Exercise 4. Make a hard- 
copy-listing of your program. This should be turned in to your . 
instructor. 



3. THE TRAPEZOID ROLE 

■ ■ ■ . '-^ ■. ., ■■ ' 

Que way to see why the left-rectangle method is some- 
what, cri^de is that tW' approxim^ltion- on the j^^ subinterval 
make^ use of only the; \laiue of f(x) at the left endpoint • 
Xj J, Although this is not too bad if the graph of f(x) 
'looks like- the one shown in Figure J, 'it is not all that 
;gp</d if thV^^i'aph .of f(x) looks like the function in 
P^^are 3. / ' ^ 3 



Figure 3. Figure 4.. * 

One way in which the error might be reduced is 
illustrated in Figure 4. Here^^we approximate 

by the integral of the function whose graph is the 
straight line connecting the points (x._^, f(x ^)) and 
(x j , f(Xj)). As Figure 4 shows, this would mean approxi- 
mating th^ area under the grapR'of f(xl over the subinter- 
Y^^ by the area pf the shaded trapezoid. This 

is the basis for the m^ethcfd k|iown -as the trapezoid rule. 
Actuairy, this and other methods of numerical integration 
that- we shall discuss follow the same approach as the' 
left ^nd' right rectangle approximations by trying to 
approximate the area under the curve of f(x) for each 
subinterval. The idea behiiji the trapezoid -rule is to 
sum areas of trapezoids ii)^f^d of arJas of rectangles. 

In general, a\trap^zoid "an be formed by the 4 points 
(Xj.p 0),^(Xj^^, f(Xj,p), (Xj, 0), (Xj^'f(Xj)) as 
illustrated by the following diagram for t\]e j^^ subint'er-* 
val (Figure 5)'. (In Figure 5 'we are really using a linear 
function, 'namely the one represented by the line through 
^C^j-l) and f(x.) to approximate the curve of f(x) in the 



.th . 

j subinterval.) From elementary geometry, we know 
the area of the trapezoid is given bv* 



A = 1/2 [V(x.;^), + f(xJ)Axt 




^ Figure 5. 

3 

Now suppose we have partitioned the interval 

a < X < b into N equal subintervals of length Ax = (b-a)/N 

as before. Thea^ we will have N trapezoids .s imi lar to the 

one described in Figure 5. The total ^area of these^N 
trapezoids' over 



■N 



a = Xq, (a + Ax) ^ x ^ , . , b = x 
may be expressed as ' • 

Ax{l/2(f(a3 + f(a,+ Ax)) + l/24f(a ^Ax) + f(a + 2Ax) ) 

(3) ^ l/2[f(a + (N-l)Ax + f(b)) } 

= Ax{f(a + Ax) + f(a + 2Ak) + ... + f(a + (N-^Mx)} 

. + i/2(f(a)' + f(b))Ax. ^ * * 

Let us find an approximation for /q1/(1 + dx by 
means of the trapezoid rule (5). Again we use N = 4 so 
Ax = 1/M and the partition is 'a =*0, x^ = 1/4, = 1/2, 
X3 = 3/4, b = 1 a«^whqin we used, the left- rectangle' me thod • 
for this example^ The approximation (3) gives 
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l/4af(}/4 * f(l/2) + f(3/4)] + 1/4 -.l/ZIfCO) + f(l)] 

= 1/4116/17 +,4/5.+ 16/25] + 1/8(1 + 1/2] 

f 4/17 +.1/5 + 4/25 + 3/a6 

= *.2353 .2000.+ .1600 +' .1875 

• = .7828. * ' 

Compare thij.s resuK for thds saJne problem with the left- 
T^ctangle method and its result as shown on page 5. 



Kxercise 5 . • , . 

Compute /j 1/x (dx) by hand calculations using the trapezoid 
rule' (3) where N » 6. Note th§ exact value In (2). 

Instead 6f w.riting a program to' compute the approxi- 
matipn by Formula (5), we shall try to write" this in an- 
other form. You should recognize that (3) is just the 
right-rectangle approximation 

(2) - f(b)Ax + lA2(f<a) f(b)]Ax. 
Now ^he quantity 

f(b)Ax + l/2(f(a) + f(b)]Ax} 
can be simplified to^ 

l/2(f(a) - f(b)]Ax. • " . . • ' ^ 

This meanj^ that if the approximating sum (2) computed by y 
the right-rectangle rule i3 adjusted by addition of the 
term^ 1/ 2 (f (a) - f(b)]Ax then we have the trapezoid rule. 

Thu^ the trapezoid rule is alJout as easy to program . 
as the rectangle rules. In fact we can -talce the program 
for computing the approximating sum by the' right-rectangle 
rule ^nd insert additional instructions*- that compute 
l/2(f(k) - f(b)]AX and then add this term to the right- 
rectangle approxiioat ing -sum. A conjjflete program for 
doing this' is now given. Does your solution, to Exercise 
2 resemble this? • ^ 

11 



100 DEF FNA(X) = 4*Xl^3 - 3*Xn + 2*X - 7 
110 INPUT A, B, N 

120 LET D = ^B-A)/N , ' ' ' j 

130 LET S = 0 

140 FOR X = (A + D TO) TO B STEP D . ^ * 

150 LET S = S + FNA (X) 

160 NEXT X. 

170 LET SI = S*D ' 

175 LET 52 = SI + D*(FNA(A) - FNA(B)-)/2 
180 PRINT S2;-"IS TRARE'XOID APPROX." 
185 . PRINT Sl; "IS RIGHT RECTANGLE APPROX." 
190 END 

Here we use line 185 to print out SI, the right- rectangle 
approximation so we fan compare this with the result from 
the trapezoid method. 

Exercise 4 , 

Type in Example Program 2 at a computer terminal. RUN this 
program for fixed values of A, B and Increasingly laraer values of 

^N. Copy down on notebook paper t*ie results from .these runs including 
your input values of A, b\ N,,/What is the smallest yalue of N that 
gives the exaot value of the integral for A B =^1? How does ' 

the trapezoid rule compare with the previous methods? Turn in your 

•solution for this exercise" to your instructor. 



Let*s review what .we have done so far in computing 
approximations to definite integrals. Our .first method 
was derived directly from the definition of the definite 
integral , The left and right rectangle approximations 
consist in approximating our given function for each 
small subinterval by a'.constant (eithej: the value of the 
function at the left eridpoint or its value at the right 
endpoint). The approximation to f(x) obtained by, piecing 
together all of these constant .functions ior a gitven 

L2 



subdivision of the interval^ la, b] is called a step 
function, and the corresponding approximating sum (left 
or right) is actually the integral of this step function. 
TKis is illustrated by the following diagram with a .sub- 
division of 5 points (Figure 6). 




Figure. , ' 

jThe trapezoidal rule us'ed a linear approximation to 
f(x) on eachr^sabinterval . The approximating function 
obtained by piecing together all of the linear pieces for 
a given subdivision of (a,b] is called^ piecewise linear, 
-and the approximating trapezoidal sura' is actually the 
integral from a to b of such a pie^ewise linear function. * 
The previous exercises should have demonstrated to you 
that the trapezoid method gives a better approximation for 
each N than the rectangle methods. 

This illustrates how complicated functions are ■ 
approximated ^n suitably small portions of their, domains 
by simpler functions, and the desired anily^is is done 
with these simpler functions. Polynomial fiinctioias are 
the best examples, of simp*le functions bec'a'us^ their values 
are easily computed and. there are simpler foriftulas for 
integrating and differentiating these. Actually, a- 
constant function is a polynomial function^of degree o 
and linear function is a polynomial function of degree one; 
thus, we might expect greater accuracy (or better fit) by 
using second-degree polynomial functions, such as parabolas. 



to approximate f(x) on suitably small subintervals and 
piece these together- to approximate the integ-ral of f(x) 
over (a ,b] . , ' 

^ ' ' * , 

4, ^SIMP50N'S .RULE ' 

' * — 

This approach leads to our ne^t method of numera.cal 
integration which is called Simpson'.s Rule and approxi- 
mates the area under the curv^ of f (x*)' over tVo adjacent 
subintervals by means of a parabola. The" approximation 
for the integral . * 

fX 



■1 



^j-1 



f (x) dx 



is derived b^ passing a parabola §j (x") through the three 

points (Xj_^,'f(x._p)., (X., f(x.)),^ andMx^,^, f(Xj,p) 

of the adjacent subintervals (x •i,rx.] and (x., x. ,] 

3^1' 3 J j+1 

The approximation to the area under fhe curve. of f(x) over 

the' adjacent subintervals (xj_^; x ^ ] and (x^ , x.^^] is 

given l?y calculating the area under the parabola g . (x) 

over tJiese subintervals. This is shown in Fi^ure\7. 




x.^^=x.+Ax 



Figure 7. 
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,0f course the overall approximation /^f (x) dx 
consists o-f the sum of the areas for these"* approximations 
by parabolas whereo j funs through the indices 1, 3, 5, ' 

N-1. A pair of adjacent subintervals' must be used 
as 3 points ^are required to' uniquely determine each approxi 
jnati-ng parabola g ^ (x)' and' there ifiust be an even number of 
> subintervals altogether so these earn be paired off. This 
meafis choices |or the yaluQs of N which represents the 
total number of sublntervals must .be even. To d'eriV.e Z 
formula for|impson's Rule to approximate /^f(x)dx*, we 
Vi\ll assume there are an even "number of subintervals^ in^ 
the partition of* (a,b3 and that each of these sub.i'ntervals 
has the same length which we den'ote, 1)y Ax as in previous ' 
discussions. , Observe as" Figure 1 indicates that x: ^ = 
- Ax and X. = x. + Ax. 

First, let us'derive a formula for the approximation 



.to 

- - +Ax 

^ f (x) dx 
'Xj -Ax 



* V" 



using these ideas. Let us as^uifte the parabola. g . (x) has ^ 
the form»gj(x) = ax^ + bx + c. The area under the pa^ra' 
bola is . ' V ^ y 



f 



x.+Ax 2 

^ (ax + bx + c)dx 
X -Ax 



1/3' ax^' + bxV2- + cx 



Xj+Ax 
Xj-- Ax ' 



Sy substitution an^d algebraic simplification, this equals 

(4-4) > l/3.Ax^(6ax.^ + 6bx. + 2aAx^ t 6c). - ' 

J J ' > ^ 

A computational* formula should involve the values 
of the original fuilttion r(x) at the point's ^x . Ax , x"*-, 
and Xj + Ax; this will occur through the determination *of 
j:he coefficients a, b;. and c. Finding the solutioy of the 
three equa^Mons ^f(Xj_j^) = g(Xj^^), f(xj^.) gCx^'il^ind 
^f^j+P"" Sf^T + l) is one technique for calculating'^.^ b. 



, and c. Appendix con-tains ' a more detailed discussion of 
this ^nd the general- method'of determining coefficients 
of the polynomial p(x) = a^^x" + a^^^^x""^ a^x + a. 

which approximates fCx) and passes through n+1 points 
(Xj, f(Xj)), (Xj^^j, f(x^^j;)) of the curve for f(x). 

^ « * 

Kn easier approach is to observe that 

g^<Xj - A^) a(Xj - Ax)^ - b(x^ - Ax),+ c, 

= a(x^ - 2Ax Xj + Ax") + b(x - Ax) + c 
and likewise — ' ^ 

'\ ' ,^j^^J ^^^^ " ^^'^j^ " ^^^^j ^^^^ M'Xj ^x) + c. 
The"re'fore, - 

gj (Xj - Ax) + gj (Xj + Ax) 

= a(2x^^ + 2Ax^) + b(2Xj) + 2c 

so tha»t if we add 

.4g. (X.) = 4ax.^ + 4bx- + 4c 
• J J J 

'to this,^ we obtain the expression' enclosed, by parentheses 
in Formula (4-a), i.e.. 



• gj(^j - ^x) * 4g. (X.) + g. (xj +*Ax) 

This means 

,^ ' - -rx.+Ax 
% -J ^ g.(x) dx 

^ -'Xj-Ax J 



2 7 ^ 

6aXj + 6bx: +2aAx + 6c. 



1/3AX (g.(x. - Ax) + 4g.(xO 



+ gj(Xj + Ax). 



Finally, 

rX.+Ax ^. 
t^-D) . J g (x) di 

''X . -AX 



■X . +Ax 

\ 

AX 

= 1/3AX (fCXj - Ax) '+ 4f(Xj) 
+ f (Xj + ^x)) . 

21 - , 
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Since^g.(x. - ^x) « f(x.'-^,Ax), g.(x.) = f(x ), and , ^ 
gjCXj + Ax) =vf(Xj + Ax) by 'construction. . ' 

If we.ad'd up the areas under these approximating 
parabolas for all rionover lapping pairs of .adjacent sub- 
. intervals,: i.e., [x^, x^) with Ix^, \] . tx^, x.) with , 

;^IX3' ^4'' 1^4' ^5^^^^^ t^5' ^6l' '-^ f-\-2'"^n-ll ^^^^ 

'^.^n-1* ^n'* ^® obtain the sum of these approximating 

,.careas from (4-b) as 

^ Ax/3[f (Xq) + 4f(x^) + fCx^)) 

^ Ax/3Tf(x2) + 4f(x ) + f(x.)) 
(5) r . ^ 

* + . . . 

' »■ » 

* An alternative way to write Formula (SJ is 

fb ' • 
f(x)dx =-Ax/3[^f(xQ) V 4f(x^-)-'-if4^-'-4-f^-^ 

• ■ * * ^ 2f(x^.2<^ 4f(x^.^) . f(x^))., 

Let us apply Formula. (5-) ' for Simpsbn's Rule to 
' at)proximate f^l/a * x^.)dx. For comparison with the left- ' 
rectangle and trapezoid methods, we" will again use N = 4 
so Ax = 1/4 and the partition is a = 0, = 1/4, = 1/2, 
Xj = 3/4, b = 1 as before. Now Formlila. (5) gives 

r/4 • l/3(f(0) + 4-f(l/4) * f(l/2)] 
>b 1/4 •. l/3[f(l/2) +^4-f(V4) * f(l)] ' ■ 

= 1/4 • 1/3[1 + 4' •- 16/17 + 4/5] ' 
*.l/4 • l/3[4/5 + 4 • 16/25 + 1/2] 
. * = l/3[l/4 ^ 16/17 ^ 1/5 >'V5 ^ 16/25 ^ 1/8) 

= l/3[. 2500 + .9412 + .'2000 -f .2000 + .6400 + .1250) 
= 1/3 [€'.3562] 

This is illustrated by Figure 8. 
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Figure 8* • 
• • t 

Figure 8 indicates how closely the parabolas' g (x) 
and g^Cx) approximate the curve y = l/(l+x^) over the 
two pairs of subinterval> and the computations illustrate 
that for about the same amount* of^effort as the trapezoid 
formula, Simpson's Rule gives much better accuracy! 
Another example which compares Simpson's Rule and the 
trapezoid .rule is given in Shenl<, p. 375. 

Consider now a function whose gra^h over the pair of 
adjacent subihtervals x.] andjx x. ^] behaves 

^as /ndicated by Figure 9(a). Here :th9 graph of the 
'quaVratic function p(x) through, the points corresponding 
'^°'''j-l' ''j' ^j+i is a poor ^pp-roximation to the 

graph, of f(x),. This. is an excel^ie_nt illustration of how 

can be made incVeasingly Accurate by 
successive doubling of the number of points^'in th^^.^ 
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y = p(x) 



y= f(x) 



y= f(x) 




Figure 9. 
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subdivision of the interval (a,b] of integration (this 
procedure* wil 1 be discussed in more detail in the next 
sectioo). As Figure 9(b) illustrates, when the original 
pair of Submtervals is^ subdi v ided intg 2 pairs of sub- 
intervals, then^a pair p^(x), P2(x) of quadrai;1^Gs functions 
is used to approximate f(x) and gives a closer fit to the 
graph. Note that* in this example, the trapezoid method 
would still not give as good an approximation, this can 
be seen by'drawing four line segments to connect pairs 
of adjacent* points on the graph of f(x) since as you 
should recall, the trapezoid method would use such 9 
piecewise linear approximat ioIrT? 

Now to implement Formula (5) bn a computer, we\ieed 
a computational routine whereby x ranges through the^ ^ 
njiddle point? of each pair of nonoverlapp mg adjacent sub.- 
ihtervals so that thb quantities f(x - Ax) + 4f(x) + ' 
f(x + Ax) can be computed and summed. To be precise, x 
must range through the values x^, x^, x^, ^ of 

our partition since tjie pairs of subihtervals are (x^., x, ] 



with- (x. 



x^] 



(x 



with. tx^_^, Xj^) 
a partition means x 



2, X3] with 1x3, x^), (x^^2% ^N-1^ 

Now using the uniform width Ax for such 



= a 



Ax , 



^N-1 



+ (N-l)Ax, i.e. 
+ (N-l)Ax in steps of 2-Ax. 

2) - 



^^3 = 



a + 3Ax, X( 



; = a + 5Ax, 
X varies from a + Ax to 
Then we can use a summing 
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vairable S to add the quantities f(x Ax) +''4f(x) +. 
f'(x + Ax) for^ach* x. Finally, according to Formula (5)',. 
S is multi^ylied by Ax/3 to complete thj computation o*f 
the sum of the areas under the parabolas. 

. Obviously, to have a computer routine for this 
calculation process would u.se a F-OR-NI:XT loop. Again 
using D for Ax, we 'have the following Rxample Program 3" 
to calculate the approximation to the* integral /y{^FNA (x) d^ 
using Formula (5) for Simpson's Ruie,^ Mere again FNA(x) " 
4x -^3x*' + 2x -'7,,b\it any function FNA(x) can be defined 
for any function by a suitab-le program instruction at line 



100. 



EXAMPLE PROGRAM o 



90 SET DIGITS 10* 

4i00 DEP FNA(x) = 4*xt3 - 3*xt2 + 2*x - 7 ' 

110 INPUT A, B, N " 

120 LET D = tB-A)/N 

1^ LET S = 0 

140- FOR X = (A + D) TO + (N-1)*D) STEP 2*D 
150 LET^S ^ S + FNA(X - D) + 4*FNA(X) 



160 >JEXT X 

( 

170 PRINT S*D/3, "IS -SlfiPSON ' S APPROXIMATION" 
18X) END ' ' 



+ FNA(X ^ 



Note the 'instruction in line 90 SET DIGITS 10 specified' 

that when the computer pri:nts njiimerical output thiis will "'^^^^ 

consist of'" 10 di§it numbers.' ^ - ? ^ 

-\ ^ '? 

Exercise 7. 



Type in Example Program 3 at a computer terminal,' Run this 

' - " - fi> 
program for fixed values of A, B and values of M like 2 4 8 16 

* ' a - ' * ' ' ' 

Remember N must be even for Simpson's Rule to worlc. Recall that the 
additional instruction given as 175 ' VRIHT BU^ -^i<y + '^2 - r*B 
- Af4 + Af3 - Af2 +^ 7*A will print the exact yaliie ;of th^. integraf ' 
/^(4x^ - 3x^ > 2x - Z)dx. 1 - \ 



Copy dovm on notebook paper the results froin your computer 
run including your, input values of A, B,, N. tOiat is the smallest 
value of H that gives the exact value of this integral for A - 0, 
B 5= 1? How does this compare with the trapezoid rule?^ Turn, in you^ 
results for this exercise to your instructor. 



An Application 

We now give an application of these approximation . 
techniques to a practical problem. Recall that if an 
object is moving at velocity v(t) as a function of time 
t, then the distance traveled between times t*= a*and 
t = b is •given by ^ 

rb 



fb 



f(t) dt , - , 

provided v( t)^^. ^ • 

"For example, silppose a motorist on a two hour trip 
nt>ted his speed at 10 minute intervals as 0, 57 , 51 , 55 , 
0, 62, 60, 53, 35, 60, 33, 35, and-O miles/hour. We*can 
es^timate the distance traveled by using either, the 
Trapezoid method or Simpson *s method (since there are an 
odd number of e.qually spaced speed values) to approxi- 
mate " . • ^ 

. l^^M dt ; ' 

here' t is time measured in hours and 10 minutes*-^ ^At 
"= 1/6 hour. ' . . ' • . 

V ' ^ 

If w6 uSe Simpson's method, then by Formula 5 we 
have - • . . 

1 ^ • * ' i' 

jf^ (0 + 4.57 > 51) ^ ^ (51 + 4-55 + 0) 

^ 3T5; (0 ^ 4- 62 + 60) + (60 + 4-58 + ^5) 
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^ 3^ (^^! 



j * 

~ 93 miles. 



+ 4-60 + 33) ^- ^ (33 + 4.35 + 0.) 



5. DOUBLING THE NUMBER OP SUBINTERVALS 
IN PARTITIONS OV UNIFORM WIDTH ^ 

,So far W6 have discussed three tecjiniques of 

numerical int^egration an^ applied each of" these to 

3 2 " - , * 

F(^c) = 4x - 3x + 2x '7 over the interval (0,1] to obtain 

computer calculated approximations. Beca-use we 'could find 

• the antiderivative for this^ example , we could compute the 

exact value of the integral and compare how close the 

approximate values were when we used different numerical 

^ techniques aiid larger value's of N. .Actually when an' anti- 

deryiative can be determined, there is jreally no point in 

using* techniques of numerical integration. On. the oth^r 

hand, sujipose we want to evaluate 

f2 ^.x2/2 • - -m ... 

, • ''2 

for which an antider ivative can.ng^^be calculated. Here 
we can not find the exact value of this integral by using 
antideriyatives.. . - - - " 

We Kave seen that larger values of N usually give 
better approximations by the numeri'caO. .methods and we 
expect this to be. true theoretically since we aro using 
Riemann Sums where x approaches 0 as N approaches infinity 
(see Shenk, p. -207). In actuality though, a machine 
cre,ated error called "round-ofr' error will build up as N « 
i-ncreases and more arithmetic computations are performed ' 

• by the computer since each number is represented 'by a 
finite number of digits in the .^comput.er and arithmetic 

'computations are 'rouhded off xiuring' ma.chine calculations-. 



27. 
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If we neglect the effect of round-off error which 
wilL be insignificant on most computers when accuracy 
to only 5 or 6 decimal places is- desired, trhen we .can ^. 
simply double the value of N each time. Thus, we would 
compute » sums $2 , S^, Sg , S^^, and when the 

difference betw&en the value, S^, of the approximating 
sum for some N and the value, S^j^, of the approximating 
sum for 2N becomes less than a specified tolerance, then 
the computational process woujd be terminated. 

Suppose for example we want an approximation accurate 
• to 0.00005 for the. integral of 

over 1-2,2]. Here we use Example Programs 2 and 3 with 
DEF FNA() instruction changedi* to DEF FNA(X) = EXP(Xf2/2) 
and Wake use of the built-in function EXP() for e^-^. 
The following tables give the results of some computer 



runs . 
N 


Trapezoid Approx. 


" N 


Simpson's Approx, 


2 


2.270671 


2 


.2.847114 


4 


2.348397 


4 


2.374305 


8 


2.581347- 


8 


2)392331 ' 


16 


2.389759 




2.3922331 


32 


2.391871 


32 


2,3925*72 


64 


2.392399 


t- 64 


^ 2.392575- 


128 


. 2.392531 






256 


^2.392565 ^ 






512 


2.392573 , 







We want |S2j^, - Sj^| < 0.00005 so 2N = 32 works for 
Simpson's Method while 2U = 256 is necessary for the 
trapezoid methpd. Here, we observe that 2.3925 is the 
approaimation for 



i. 



'-2 " , 

that "is accurate to O^^OOOOS, i.e., 4 decimal places,- $.nd 
is obtained by both methods. 
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Observe that the example programs can be used lo 
compute approximations for the integrals of other func- 
tions by simpiy changing the DEF FNA() instruction. For 

'F(x) = e"^^^^ 
we used tHilr built- in function EXP(); other built-in 
functions are SIN-() , LOG () -natural logarithm, COS()^ etc. 
Appendix A contains a "complete list of such built-in 
function's available m most versions of thp programming 
language BASIC. \ 

; / 

Exercise 8 

Use th^ trapezoid and Simpson techniques to compute sin(ln(x)) dx 
accurate to 0.0000005 by doublin*g N for successive computer runs. 
Here DH^ FNA(x) = SIN (LOG(x)). What are the val&es of N required 
for each technique? Remember to insert the instruction 5 SET DIGITS 
10. 

Exercise 9 . Since L^l/(l+x^)dx - Arctan (1) = 77/4 then 
fl 2 

4 • /qI/CI + X ) dx =77. Use the method of doubling the number of 
subinterv^ls to find the value of ir accurate to 8 decimal places. 
What are the values of N rec(uired for this technique? Use only 
Simpson's Method (Example Pifo'^gTam 3) and the additional instruction 
175 PRINT 4*S*D/3, »'IS THE APPROXIMATION TO ir". 
Note: TT = 3.14159265^58979323846264... 



In practice, most numerical integration techniques 
are carried out on a computer; and the accuracy of the 
approximation is usually determined by doubling the number 
of subintervals until two successive approximations differ 
by less than :the prescribed tolerance. This "s topping^!!^ 
rule", although widely used, does not always work 
efficiently. Consider the Figure 10; this shows that 
unless the number of subdivision points is large enough 
to^ require the integrand to be evaluated at points in the 
intervaT [c, c+{i] very little change will^be noticed 
betvj^een two successive approximations. 



00 
1^ o 
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"T 1 1 1 • 

a c C+h b 



Figure 10. 

There is another reason for doubling N until two 
successive approximat ions^di f f er by less than the pre- 
scribed tole^rance. If an "optimal'' N were on the order 
of 40*0 for example, one would expect S^qq and S^^^ to, be 
• v^T'y close to each other, even if neither were close 
enough to the correct approximation. Actually, a more 
elaborate computational Algorithm can incorporate the 

computations used in the. calculation of S,, when S^., is 

N . 2N 

computed. This approach would give a more ^ef f icient 
program and eliminate unnecessary duplication .of compu- ' 
tations which occurs when a straightforward program like 
Example Program 3' i's simply repeated for N = 2 , 4,8,.... 
An example of such an algorithm for Simpson's Method is 
described by D.A. Smith (21, ^p. 121-122). 

It is important to note that Simpson/s rule requires 
an evenly spaced parti t ion with an even number of subin- 
t-erv^ls. Thus, given y = F(x) if the x- values do not 
satisfy this criteria then the trapezoid method is the 
only one of these two methods' that can be used. -This 
will only occur for . functions defined from experimental 
data as in' the appl icat ions example of Section 4. When 
^.^-^^V'^.^VP'^ is def ined by a_ fo_rmula,,^,.then any given inter- 
val of integration can be partitioned in such a way that ' 
Simpson's rule can be applied. Several ^xamples have 
demonstrated that foV N = 8, accuracy to 4 decimal places 
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can usually be obtained by Simpson's method.. This 
means that hand computations could suffice if a computer 
is not available. Because of the preceeding observa- 
tions, it should not be surprising that Simpson's rule 
is widely used in practical computations. 
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. ■ -t-- "^ ^ ^^ 'The 'i&Ji^^^^^ ' V 

functions,' tWJ caICl^^^^■?^^i^^^ of aa^,- : ^ - r 

^'-eprap laxity and- laiiy" ilficf u^^^ Tbe ^ 
quantity; CneX i^' aiv arg^^Jit^-tpaj^iH^^ 



function :;. ' D£CSfel PTlp>^-v- ' '/ / ^- ' ^ ' ' 
AES(ne-)-:L. >..r/\'rJEfnda- -ther^ab Solute v^l 5b ^ u^.^'v'' 





ATN(ne:D 


,"3Sjid"S''tKe.:'at^in^^ rti the principal^ 
value ^ng.e ' VttA 't'd— +tt72 > - ' -v- 


• • 


COS^Cne) . S' 


P in d s ^ )ve *c o ai n e- ^c5^1 jiie5?''::the' "dri^^fi"^^' -is 
-expirees s e J'^in^*. x.^^i^tj^it^v:'^^;; / - - . --^ "r^ 




EXP(ne) • 


Fiad^'-'tliiT^vaJiie-"':© to ftHe' pbHtU". pJi'j^p , 




l.NT(»ne) 


■ finds t he '^aT-g'e&t-. -integer- not greater t^an' 
ne. - -^Ex^jnprer^ IN.T(5.95J '=''5^ and IN.T(-5.9S) 
= -6. - * ' - -.r - \ - 


\ 


^LGT(ne) 


FjLnds the base io 'XbgaritTim of ne; ne>0jjj- 
otherwise,.5Ln: execution error causes program 
termination, ' ' . 




• LOG(ne) 


Finds the na.turaT logarithm of ne; ne>0, ^ 
otherwise an execution error causes program 
termination. \' * . * ' ' " . 




^ SIN(ne)^. 


Finds the sine* of .tie; the' angle ne is ex-r 
pressed-in radians. 




SQR(ne) 
*^ « * 


Finds the. square root of ne; ne > 0, other- 
wise an executijon error causes program 
termination. ' . ^ 




TAN(ne) 


Finds' the tangent of ne ; -the 'angle ne is 
expr.essed in radians. - 
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APPENDI-X R 



:5upp^se we want to -app'55^4mate the cdrve* of ,iF(x) - 
"T/l^M by a quadratic .function'g fx) - ax + bx + c qvet 

the. interval [0,1].. " must/il^b three points xSn t he 
> graph o^(x> say: (1,475), and p, 1/2). These - - 
c hoi c e^^na k e " US o f - 1 h'^^-- p ar t i t i o 1 0 , 1 / 2 ) , ( 1 /i^ T p q f ' ' 
;the.-inteFvai [0 , llV r Si^ice g'<6r=;F(0), gn/2) ^.F(1/2T, 
'rind gClJ '= F(l) we .have- the s.)t^t^ t^f 3 equation^-vi th 
uhicndwns^ a^,- b-, c-^s follaws:V ^ . 

a^..oi.+ b-0 1_ * / . ^ c = 1 



^aClv^2) '% b(l*/r) + c f 4/5 a-1/4 +-b-V2 + c = 4/5 
.a'(lj^ + b(l) '- 1:^-,= 
Th'is simplifies to' 



• • ' 2 ' ' 

a'(lj + b(l)-t:^-,= 1-/2 a.l + b-1 + c = 1/2 



c = 1 

5a + 40b + 20c = 16 

2a * 2b + 2c = 1 • 

The solutions are c = 1, a = -1/5, and b = -3/10 so that 

. " g(x) = - 3^ . 1.- ' • 



Observe that 
rl 



g(x) dx = 

^ 



il 
'60 



and Formula (4-^) of Simpson's approximation gives 

.1 ^ ^ \ - 

J^F(x)dx- i . i(F(0) + 4F(i/2) >.F(1)) = ^ 

also when the subdivision 10-= x^., 1/2 = x^, 1 = of 
[0,1] is used. ^ ' ' 

The above procedure is often called the method of 

undetermined coefficients. If we had 4 points of some 

fimcfion, iay (:J,l/2), (6,1), (1,2), (2,4)', we could 

Or: ■ 



4 ' ' 

approximate this by a third degree' polynomial p(x) 
ax^ + bx^ + cx + d: The equations pC-1) = 1/2, ' 
p(0) = 1, p(l) = 2, and p(2) = 4 translate into 



+ * b c + d = 1/2 



d = 1 



a+b+c + d = 7. , ' 

8a + 4b + 2c + d = 4. 

After some labor, one finds a = 1/12, b'= 1/4, <; = 2/3', 
d = 1, and • 

\ in general," to approximate a function by a poly- 
nomial of degree n, i.e., ^ 

p(x) = a^x^ ^ a„_,x""^ + a 



n-l^- ^ ... - a^x > ?Q 



requires (n+i) points on the graph of Jp(x) , These will 
give n+1. linear equations involving the n+1 unknowns a , 
^n-1' ^n-2''*''* ^0' above examples illustrate, 

the method of undetermined coefficients can become quite 
tedious. A simpler method called Lagrange Interpolation 
has berfn devised. A good discussion of this may be found 
in Flanders (5, p,377>379). 



STUDENT FORM 1 
Request for Help 



Return to: 

EDC/UMAP 

55* Chapel St. 

Newjton, MA 02160 



Student; If you have trouble with a specific part of this unit, please fill 
out this form and take it to your instructor for assistance. The information 
you give will help the author^ .to revise the unit. 



Your Name 




\ 


Page ^ 




r 




O Upper 


OR 


' Section 


OR 


QMiddle- 




Paragraph 




Lower 









Unit No, 



Descripti<^n of Difficulty: (Please be specific) 



.Model Exam 
Problem No, 

Text, 
Problem No, 




Instructor t Please indicate your resolution of the difficulty in this box. 



o 



Corrected errors in materials, y^t^t co^ections here: 



o 



Gave student better explanation, example, or procedure than in unit. 
Give brief outline of your addition here: 



o 



\. . ■ 

Assisted student in acquiring general learning and probl^-solving 
skills (not using examp^^s from this unit.) 
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Instructor's Signature_ 



^Please use reverse if necessary. 



Return to: 

STUDENT FORM 2. EDC/UMAP 

^ Ti 4. A ,.4 4 ' 55 Chapel St. 

Unit Questionnaire ^^^^^^^^ ^ Q2160 



Name Unit No. Date 



Institution Course No. 



Check the choice for each question that comes closest to your personal opinion. 
1. Hov -useful was the amount of detail^^n the unit ? 
Not enough detail to understand the pnit 

Unit would have* been clearer; with more detail \^ 

^Appropriate amount of* detail 



JJnit was occasionally too^ detailed, but thi^ was not distracting 
Too much details ly^as often distracted ^ 



How helpful were the problem answers ? 

S ample solutions were too brief; I could ^tf^ do the intermediate steps 
^Sufficient information -was gi^fiH to solve » the problems 



^ S^ple solutions ^were too detailed; I didn't need them 

^ . . 'l"N 

3. Except ^or fulfilling the prerequisites, how^much did you use other sources (for 
example, instructor, friends, ^ other books) in order to understand the unit? 

A Lot S omewhat A Little Not at all 

A. How long was this unit in comparison to the amount of time you generally sp^d^on 
, a- lesson (lectute and hom^crfk assignment) in. a typical math or science course? 

Much Somewhat About , ' Somewhat Much 
* ' Longer Longer • , ' t he- Same ^Shorter ^Shorter • 



5. Were" any of the following »parts of the unit confusing or distracting? (Check 
as many as apply .),"/' 

Prerequisites ^ ^ - 

Statement of skills and concepts (objectives) 



^Paragraph headings 

^Examples ^ ' . 
[[Special Assistance Supplement (if present)^ 
jOther, please explain : 



6. Were any of the following parts of the .unit particularly helpful? (Check as many 
as apply.) - . . 

' P rerequigites 

^Statemejit of skills' and concepts (objectives) 

" E xamples «- ^ . ^ . ^ 

Problems ' . • ' ' ^ ^ 



__P»ragraph headings 
Jlable of Contents ^ 

^Special Assistance Supplement (if presentX 
"other, please explain i 



Please describe anything in t>he unit that you did not particularly like. 

Please describe anything that you found particularly helpful. . (Please use the back pf 
this sheet'~if you need more space.) \^ ^* 



